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1. Introduction 

The only way to remove pending issues associated with systematic errors due to quenching 
in lattice QCD is carrying out simulation with dynamical quarks. Among several ongoing projects 
generating 2+1 flavour dynamical gauge ensembles, a striking characteristics of that by the RBC- 
UKQCD Collaborations is in their usage of domain wall fermions (DWF) for the sea quarks. DWF 
breaks chiral symmetry only by a small amount and thus has least discretization error. There are a 
lot of extensive measurements planned on the ensemble in near future. 

We report our calculation of the static quark potential and preliminary results of lattice spacing 
a from Sommer scale ro ~ 0.5 fm on the RBC-UKQCD gauge ensembles. 



2. Calculation 

The static potential V{f) between a quark and anti-quark pair at relative spatial displacement 
r, is calculated from the Wilson loop {W{r,t)): {W{r,t)) = C(r)e"^('''' + ("excited states"), with a 
normalisation C(r = 0) = 1. 

We calculate (W^(r,f)) in lattice QCD with 2+1 flavour domain wall fermion |Q] and several 
gauge actions (Iwasaki [^], DBW2 iQ], and more negative rectangle coefficient, ci, actions) on 
16^ X 32 lattice with fifth dimension size, Lv=8, and domain wall height, M5=1.8. The reason for 
calculating several gauge actions was to search for actions producing smaller residual mass, 01^5, at 
fixed Ls and lattice spacing a [^. In quenched simulations, changing c\ caused a drastic reduction 
of /Mres: wii-^''**""' 10 X mf^^^'^^^ 100 X mr™*^^ at w 2 GeV [^]. We use three degenerate 
up and down quark masses, mu^d = 0.01, 0.02, 0.04, and the strange quark mass = 0.04 was our 
best pre-simulation estimate of the physical strange quark mass. Simulation parameters are tabled 
in Table [I[ I and |. 

We compare two independent analyses, both of which implement APE smearing [Q] for spatial 
links to improve the signal/noise. 

The main difference between these analyses is the path connecting the quark and anti-quark 
in a time slice (Figure [T]). The first analysis (Type-I) is carried out using the closest path to the 
diagonal line connecting the quark-antiquark. This path is determined by Bresenham algorithm [^] 
in the same way as [|, 10]. The second path (Type-II) is the taxi driver's path employed by Chi^oma 



code [11]. The smearing coefficient and number of smearing steps are tuned to be (csmeai jA'^smear; 
(0.50,20 ~ 40), where the overlap with ground state, C(r), is approximately its maximum. 

The physical parameters are obtained by fitting the lattice value oiVir) to the function: Vir) = 
Vb — o;/|r| + a|r|, from which Sommer scale ra/a = (1.65 — a)/o is determined in lattice unit. 
For each gauge action and fixed gauge coupling j3, we estimate lattice scale a^^ assuming ro = 0.5 
fm in the chiral limit. Note this assumption turns out to be consistent with the other scale setting 
ansatz using rho meson mass, nip = 110 MeV in previous simulations [p, 12]. 



3. Results 

Figure ^ shows potentials measured using Type-I and Type-II paths. These two analyses were 
performed on same ensemble (j3 = 2.2, m„ ,^ = nis = 0.04 on Table |l]), and central values are con- 
sistent with each other within statistical error. Sharply impoved signals (red circles) for relatively 
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Figure 1: Sketch of the path connecting quark and " ^ 
anti-quark. Type-I (red) is the closest path to diago- 
nal line (green) by Bresenham algorithm, and Type-II 0.4 
(blue) is the taxi driver's path. 




Figure 2: y(r) vs |r| for Iwasaki action, /3 — 2.20, 
nil, d = 'w.v = 0.04. 



long distance off-axis spatial paths indicate that the Q-Q state by Type-I has larger overlap with the 
ground state than that by Type-II. This is reasonably consistent with the picture of the QCD-string 
connecting quark and anti-quark. 

The time separation of Wilson loop, t, is chosen to minimise excited states contamination to the 
ground state potential by monitoring t dependence of V(r). Then V(r) is fit to the formula for |r| G 
\>'min^T'max\■ Figure || show tq as a function of dynamical quark mass, + mies> for several fixed 
j3's. To estimate lattice spacing a, we take the extrapolation to the point »i„,f/ + wires = linearly 
and obtain ro/a in this limit (open symbols) from simulation points (filled symbols). Preliminary 
results of this extrapolation and the value of cr^ are summarised in Table |l], ^ and ||, in which 
data at the simulation points are listed as well. For the Type-I data, we estimate systematic errors 
by the dependence of the results on the fit ranges in extracting potentials from Wilson loop and 
in fitting the potentials to the specific function introduced in Section 1. We estimate these errors 
from the variation of central values for ?-,„,•„ G [\/2, \/6], rmax G [7,9], f = 5,6. Available results for 
Iwasaki j8 = 2.13 and DBW2 jS = 0.764 are obtained by Type-II path. Analysis by Type-I path is 
under progress. In particular, RBC-UKQCD has started their major production run using Iwasaki 
/3 = 2.13 on 24^ x 64 lattice with = 16. Our results for Iwasaki j3 = 2.13 provide properties of 
this ensemble. 

Finally we measure the magnitude of discretization error using the data from all ensembles 
we measured. Due to the (almost exact) chiral symmetry of DWF, the scaling violation is expected 



as ^((oAqcd) ) + &{am,, 



for the parameters used. We examine this expectation by 



plotting a pair of dimensionless quantities, {{r^ninY ,rQmp) in Figure Q niji, nip ai^e pseudo-scalai^ 
and vector meson masses respectively at relatively large valence quark mass set to be equal to 
sea quark mass m,, ^ [^]. If scaling were perfect (no discretization eiTor), all data points are 
on a universal curve in the two dimensional plane. The left panel of Figure ^ demonstrates the 
expectation is true for these quantities: points are on a line within statistical eiTor ( < 5 % ). 

We also compare DBW2 j3 = 0.80 results with those of two flavour (red square) and quenched 
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Figure 3: ro as a function of dynamical quark mass for Iwasaki (left) and DBW2 (right) gauge action. Each 
line shows linear extrapolation to the chiral limit, m„ f/ + mres — * 0, for fixed j3 's indicated in legends. 
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# conf 


Alg 


Path 




a-^ [GeV] 




0.04 


320 






4.788(35)(30) 


1.890(14)(12) 


2.20 


0.02 


177 


RHMC 


Type-I 


5.005(44)(48) 


1.975(17)(19) 




-mres = -0.0065(1) 








5.293(1 12)(81) 


2.089(44)(32) 




0.04 


306 






4.069(46) 


1.606(18) 


2.13 


0.02 


300 


RHMC 


Type-II 


4.270(50) 


1.685(20) 




-mres = -0.0104(2) 








4.576(145) 


1.806(57) 



Table 1: Results of ro for Iwasaki gauge action (ci = —0.331). m.s = 0.04, = 8, M5 = 1.8. Note our 
estimation for lattice cut off, a"', is the values at chiral limit, — > — Wres- 







# conf 


Alg 


Path 




a-^ [GeV] 


0.80 


0.04 


200 


Ralg 


Type 
-I 


5.014(38)(67) 


1.979(15)(26) 


0.78 


0.04 
0.02 

-m,^^ = -0.00437(6) 


165 
180 


RHMC 


Type 
-I 


4.600(49)(16) 
4.863(54)(142) 
5.184(134)(335) 


1.815(19)(6) 
1.929(21)(56) 
2.046(53)(132) 


0.764 


0.04 
0.02 

-m,^^ = -0.00546(7) 


310 
300 


RHMC 


Type 
-II 


4.415(41) 
4.737(105) 
5.147(244) 


1.742(16) 
1.869(41) 
2.031(96) 


0.72 


0.04 
0.02 
0.01 

-mres = -0.01066(7) 


280 
400 
150 


RHMC 


Type 
-I 


3.554(31)(67) 
3.841(39)(61) 
4.029(63)(27) 
4.319(80)(28) 


1.403(12)(27) 
1.516(15)(24) 
1.590(25)(11) 
1.704(32)(11) 



Table 2: Same as Table 1 for DBW2 action (cj = -1.4069). 
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c\ 


# conf 


Alg 


Path 


n)/a 


a-^ [GeV] 


0.53 


-2.3 


160 


Ralg 


Type-I 


5.074(48)(48) 


2.002(19)(19) 


0.48 




180 






4.002(57)(60) 


1.579(22)(24) 


0.36 




135 






5.060(66)(33) 


1.997(26)(13) 


0.33 


-3.57 


145 


Ralg 


Type-I 


4.335(47)(102) 


1.711(19)(40) 


0.32 




180 






3.968(49)(43) 


1.566(19)(17) 


0.16 


-7.47 


165 


Ralg 


Type-I 


4.050(57)(43) 


1.598(22X17) 



Table 3: Results of r^/a and corresponding lattice scale a from the simulations with various ci's with 
mu,d = = 0.04. 




Figure 4: r^nip vs. {rQm-,i)~ for Nf = 2 + \ with m^ = 0.04 (left). Right panel is the same plot for DBW2 
Nf = 0{P = lM),Nf = 2{j5=QM),andNf = 2+l (j3 ==0.72, m, = 0.04) with linear chkal extrapolation 
{r^nijif' — > for Nf = 2 and 0. Error is uncorrected between ro and meson spectra. 

(green diamonds) simulation in the right panel. The universal line of two flavour points roughly 
accommodate the 2+1 flavour points. It is expected because is heavy and has less effect on 
up and down quark meson masses than u and d quarks. On the other hand, quenched data have 
significantly smaller slope than other two data sets leading larger value of r^nip in the chiral limit. 

However, we note that the estimation of statistical errors are not so accurate that they reflect 
the actuality completely since the length of total trajectory and the separation between measured 
trajectories may not be large enough for some of ensembles. 

Acknowledgement 

We thank Peter Boyle, Mike Clai^k, Saul Cohen, Sam Li, Meifeng Lin, Chris Maynard, and 
Robert Tweedie, and Azusa Yamaguchi for generating ensembles and meson masses used in this 
work. We thank Peter Boyle, Dong Chen, Norman Christ, Mike Clark, Saul Cohen, Calin Cris- 
tian, Zhihua Dong, Alan Gara, Andrew Jackson, Balint Joo, Chulwoo Jung, Richard Kenway, 



093/5 



The static quark potential in 2+1 flavour Domain Wall QCD from QCDOC 



Koichi Hashimoto 



Changhoan Kim, Ludmila Levkova, Xiaodong Liao, Guofeng Liu, Robert Mawhinney, Shigemi 
Ohta, Konstantin Petrov, Tilo Wettig and Azusa Yamaguchi for developing with us the QCDOC ma- 
chine and its software. This development and the resulting computer equipment used in this calcula- 
tion were funded by the U.S. DOE grant DE-FG02-92ER40699, PPARC JIF grant PPA/J/S/1998/00756 
and by RIKEN. This work was supported by DOE grant DE-FG02-92ER40699 and we thank 
RIKEN, BNL, the U.S. DOE and University of Edinburgh for providing the facilities essential 
for the completion of this work. We also thank the RIKEN Super Combined Cluster (RSCC) at 
RIKEN, for the computer resources used for the static quark potential calculation. K.H. thanks 
RIKEN BNL Research Center for its hospitality where this work was performed. 

References 

[1] R. Sommer, "A New Way to Set the Energy Scale in Lattice Gauge Theories and its Application to the 
Static Force and a, in SU(2) Yang-Mills Theory", Nucl. Phys. B41 1 (1994) 839, 

[hep-lat/9310022] . 

[2] D. B. Kaplan, "A Method for Simulating Chiral Fermions on the Lattice", Phys. Lett. B288 (1992) 
342, [hep-lat/9206013] ; Y. Shamir, "Chiral Fermions from Lattice Boundaries" , 'Hud. Phys. 
B406 (1993) 90, [hep-lat/ 9303005 ] ; V. Furman, Y. Shamir, "Axial symmetries in lattice QCD 
with Kaplan fermions" , Nucl. Phys. B439 (1995) 54 [hep-lat/ 94 05004 ] . 

[3] Y. Iwasaki, "Renormalization Group Analysis of Lattice Theories and Improved Lattice Action ", 
uTHEP-117, uTHEP-118 (1983); Y. Iwasaki and T. Yoshie, "Renormalization group improved action 
for SU (3) lattice gauge theory and the string tension ", Phys. Lett. B143 (1984) 449. 

[4] T. Takaishi, "Heavy quark potential and effective actions on blocked configurations ", Phys. Rev. D54 
(1996) 1050; P. de Forcrand et al. (QCD-TARO), "Renormalization group flow ofSU(3) lattice gauge 
theory - Numerical studies in a two coupling space", Nucl. Phys. B577 (2000) 263 

[hep-lat/9911033] . 

[5] R. D. Mawhinney (RBC and UKQCD) in these proceedings. 

[6] Y. Aoki et al. (RBC), "Domain wall fermions with improved gauge actions", Phys. Rev. D69 (2004) 
074504, [hep/lat-0211023] . 

[7] M. Albanese et al. (APE), "Glueball masses and String tension in Lattice QCD", Phys. Lett. B192 
(1987) 163. 

[8] B. Bolder et al., "A High Precision Study of the QQ Potential from Wilson Loops in the Regime of 
String Breaking", Phys. Rev D63 (2001) 074504, [hep-lat / 0411006]. 

[9] Y. Aoki et al. (RBC), "Lattice QCD with two dynamical flavors of domain wall quarks", 

hep-lat/0411006. 

[10] K. Hashimoto and T. Izubuchi (RBC), "Static Q-Q Potential from Nf = 2 Dynamical Domain-Wall 
QCD", Nucl. Phys B140(Proc. Suppl.) (2005) 341, [hep-lat/04 9101 ] . 

[11] R. G. Edwards and B. Joo, "The Chroma Software System for Lattice QCD", Nucl. Phys B140 (Proc. 
Suppl.) (2005)832, [hep-lat/040 9003] . 

[12] T. Izubuchi (RBC), "Hadron Spectrum and Decay Constant from Nf = 2 Domain Wall QCD", Nucl. 
Phys B 140 (Proc. Suppl.) (2005) 237, [ hep-lat / 0410034]. 

[13] M. Lin (RBC and UKQCD) in these proceedings; C. M. Maynard (RBC and UKQCD) in these 
proceedings; R. Tweedie (RBC and UKQCD) in these proceedings. 



093/6 



